Graphene field effect transistors as room-temperature Terahertz detectors 



o 



L. Vicarelli/ M.S. Vitiello/'Q D. Coquillat,^'^ A. Lombardo,'^ A.C. 

Ferrari, "^ W. Knap,^'^ M. Polini,^ V. Pellegrini,^ and A. Tredicucci^ 

^NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, 1-56126 Pisa, Italy 

^ Universite Montpellier 2, Laboratoire Charles Coulomb UMR 5221, F-34095, Montpellier, France 

^CNRS, Laboratoire Charles Coulomb UMR 5221, F-34095, Montpellier, France 

"^Department of Engineering, Cambridge University, Cambridge, CBS OFA, UK 

The unique optoelectronic properties of graphene [l] make it an ideal platform for a variety of 
photonic applications [2], including fast photodetectors [3], transparent electrodes |4], optical mod- 
ulators [5], and ultra- fast lasers \6^. Owing to its high carrier mobility, gapless spectrum, and 
frequency- independent absorption coefficient, it has been recognized as a very promising element 
for the development of detectors and modulators operating in the Terahertz (THz) region of the 
electromagnetic spectrum (wavelengths in the hundreds of /xm range) , which is still severely lacking 
in terms of solid-state devices. Here we demonstrate efficient THz detectors based on antenna- 
coupled graphene field-effect transistors (FETs). These exploit the non- linear FET response to the 
oscillating radiation field at the gate electrode, with contributions of thermoelectric and photocon- 
ductive origin. We demonstrate room temperature (RT) operation at 0.3 THz, with noise equivalent 
power (NEP) levels < 30 nW/Hz^^^, showing that our devices are well beyond a proof-of-concept 
phase and can already be used in a realistic setting, enabling large area, fast imaging of macroscopic 
samples. 
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Photodetection of far-infrared radiation (from hun- 
dreds of GHz to a few THz) is technologically impor- 
tant for a wide range of potential applications spanning 
from medical diagnostics to process control, and home- 
land security [7 . THz radiation penetrates through nu- 
merous commonly used dielectric materials, otherwise 
opaque for visible and mid-infrared light. At the same 
time it allows spectroscopic identification of hazardous 
substances and compounds through their characteristic 
molecular fingerprints [7 . Commercially available THz 
detectors are based on thermal sensing elements being 
either very slow (10 — 400 Hz modulation frequency) 
(Golay cells, pyroelectric elements), or requiring deep 
cryogenic cooling (4 K) (superconducting hot-electron 
bolometers) ^, while those exploiting fast non-linear 
electronics (Schottky diodes) are mostly limited to the 
range below 1 THz [8]. 

A more recent approach exploits FETs [9], either 
as HTV high-electron-mobility transistors (HEMTs) 
or Si-based complementary metal-oxide-semiconductors 
(CMOSs). These devices provide excellent sensitivities 
(NEP in the 10"^° - 10"^^ W/Hz^/^ range depending on 
temperature and operation frequency), with the intrin- 
sic possibility of high-speed response (in principle just 
limited by the read-out electronics impedance). This ap- 
proach has been also demonstrated with self- assembled 
In As nanowires, allowing RT THz detection with NEPs 
of 10-^ W/Hz^/^ at 0.3 THz [10 . 

THz detection in FETs is mediated by the excitation 
of plasma waves in the transistor channel [TT]. On one 
hand, a strong resonant photoresponse is predicted in 
materials having plasma damping rates lower than both 
the frequency of the incoming radiation and the electron 
transit time in the transistor channel. This requires elec- 



tron mobilities of at least a few thousand cm^/(V x s) 
at frequencies above 1 THz. Under these conditions, 
stationary states arise from the quantization of plasma 
waves over the gate length. However, this concept was 
never demonstrated at room temperature. Furthermore, 
even though some results have been achieved at cryogenic 
temperatures [9 , thus far no clear-cut evidence was re- 
ported for the expected large responsivity enhancement. 
On the other hand, when plasma oscillations are over- 
damped, broadband THz detection is observed [9 . In 
this case, the oscillating electric field of the incoming ra- 
diation applied between source and gate electrodes pro- 
duces a modulation of both charge density and carrier 
drift velocity [11]. Carriers traveling towards the drain 
generate a continuous source-drain voltage Au, which is 
controlled by the carrier density in the channel, and can 
be then maximized by varying the gate voltage. 

A material having high carrier mobility at RT is crucial 
to take full advantage of the performance enhancement 
expected from resonant detection, and to make FET 
detectors equally performant at operating frequencies 
> 1 THz [12 . The naturally occurring two-dimensional 
electron gas in a doped graphene sheet has a very high 
mobility even at RT [1 . Additionally, it supports plasma 
waves [13 which are weakly damped [Ml |15, in high- 
quality samples. Thus, graphene FET plasma-based pho- 
todetectors could outperform other THz detection tech- 
nologies [16, 17 . 

So far, the development of graphene-based photode- 
tectors has been limited to the near-infrared region, ex- 
ploiting the creation of electron-hole pairs following light 
absorption [3, 18-20 . In the THz range, however, the 
photon energy is just a few meV, and light absorption 
is prevented by Pauli blocking from unavoidable dop- 



ing of as-produced graphene samples ^V or charge in- 
homogeneities [22 . Cryogenic graphene bolometers have 
also been recently developed and their operation demon- 
strated in the mid-IR [23^ . 

Here we demonstrate efficient, RT THz detectors, by 
combining top-gated graphene FETs with THz receiver 
architectures designed to achieve the desired device sensi- 
tivity. When a THz radiation beam is focused on a FET 
channel, the incoming optical wave is rectified, leading to 
a dc signal Au between source and drain, proportional 
to the received optical power [9 . To induce a detectable 
photovoltage signal across the FET, an asymmetry be- 
tween source and drain contacts is needed, for example 
in feeding the incoming radiation. This can be achieved 
either by an asymmetric design of source and drain con- 
tacts or by using low shunt-capacitance antennas to fun- 
nel the radiation into the strongly sub-wavelength de- 
tecting element. Antenna coupling ensures a selective re- 
sponsivity to both the spatial mode and the polarization 
of the incoming radiation. 

We use graphene exfoliated on Si/Si02 (see Methods) 
and lithographically define a single lobe of a log-periodic 
circular-toothed antenna [10 as source contact, with an 
outer radius of 322 /im, while the drain is a metal line 
running to the bonding pad. A second identical antenna 
lobe acting as top-gate contact is then defined by e-beam 
lithography, after depositing via atomic layer deposition 
a ~ 35 nm thick Hf02 insulating layer. The channel 
length is between 7 and 10 /im, while the gate length (Lg) 
is 200 nm or 300 nm depending on the device. Figure [l] 
shows a microscope image of a fabricated graphene FET, 
embedded within a schematic representation of the THz 
detection configuration. 

Figures |2Va) and [2Vc) show the conductivity (a) ex- 
tracted from the current (/sd) measured at RT in mono- 
layer (SLG) and bilayer (BLG) graphene devices, re- 
spectively, while sweeping the gate voltage (Vg) from 
— 1 V to +3.5 V and keeping the source-drain voltage 
at Vsd = 0.1 mV. In both cases, as expected [1 , vary- 
ing the gate voltage, the conductivity goes through a 
minimum when the chemical potential below the gate is 
tuned across the charge neutrality point (CNP). In both 
samples, the CNP is located at positive gate voltages, un- 
derlying the p-doping of our graphene samples consistent 
with the doping level measured by Raman spectroscopy. 
Note that in our devices the gate electrode length is much 
smaller than the source-drain distance. This implies that, 
in general, one cannot simply assume the measured resis- 
tance to be fully dominated by that of the gated region. 

A diffusive model (Ref. [24^ and Methods) of transport 
predicts a second-order non-linear response when an os- 
cillating THz field is applied between gate and source of 
our FETs. This model yields a photovoltage Au propor- 
tional to the derivative of the channel conductivity with 
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FIG. 1: (Color online) Schematics of the THz detection con- 
figuration in a Field-Effect Transistor (FET) embedding the 
optical image of the central area of a bilayer graphene-based 
FET. A log-periodic circular-toothed antenna having an outer 
radius of 322 /am is patterned between the source and gate 
contacts, while the drain is a metal line running to the bond- 
ing pad. The channel length is ?^ 10 /xm, while the gate length 
(Lg) is 200 nm. The metallic (Cr/Au) top-gate is patterned 
in the middle of the graphene channel with standard electron- 
beam lithography and lift-off technique after atomic-layer de- 
position of a 20 nm-thick Hf02 insulating layer. An elec- 
tromagnetic simulation using a finite-element method (FEM) 
software (Ansoft HFSS) to calculate the resonance frequen- 
cies of the log-periodic antenna (without the bow-tie element) 
leads to resonant frequencies close to 0.4 THz, 0.7 THz, 
1 THz, and 1.4 THz. The THz radiation at 0.3 THz is fo- 
cused on the bow-tie device by off-axis parabolic mirrors with 
an optical power of 2.1 mW and a spot size with a diameter 
of roughly 30 microns. Inset: Scanning Electron Microscope 
image of the monolayer graphene FET. 
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Figures ^b)-(d) show the predicted photovoltage as a 
function of gate bias for the SLG ^) and BLG ([2]i), 
as inferred from Eq. ([l]), using as an input the measured 
gate- voltage-dependent conductivity a = cr(Fg). In both 
cases, Au varies from negative to positive, while crossing 
the CNP, consistent with the ambipolar nature of charge 
carrier transport in graphene [1 . 

To access the photoresponse, we shine the 0.3 THz 
radiation generated by an electronic source based on fre- 
quency multipliers [25 . The radiation is collimated and 
focused by a set of two f/1 off-axis parabolic mirrors. The 
intensity is mechanically chopped at frequencies varying 



1.5x10 



^ 1.0x10 



0.5x10"' 








0.05 



Vg(V) 



Vg(V) 



FIG. 2: (Color online) (a,c) Representative electrical conduc- 
tivity plots measured at room temperature as a function of 
the gate voltage in the single-layer graphene (a) and bilayer 
graphene (c) field-effect transistor while keeping the drain bias 
at 1/sd = 0.1 mV. The drain contact was connected to a 
current amplifier converting the current into a voltage signal 
with an amplification factor of 10^ V/A; (b,d) Derivative of 
the device conductivity as a function of the gate voltage Vg, 
plotted once multiplied for the device resistance. This curve 
represents the expected device responsivity following a diffu- 
sive model of transport as described in Ref. [23] and in the 
Methods. The change of sign occurs as Vg is swept through 
the charge neutrality point. 



between 90 Hz and 1 kHz, and Au is measured by means 
of a lock-in amplifier in series with a voltage preamplifier 
having an input impedance of 10 MQ and an amplifica- 
tion factor G = 1000. The vertically polarized incoming 
radiation impinges from the free space onto the graphene 
devices mounted on a dual-in-line package with an optical 
power Pt = 2.1 mW. 

The photoresponse is measured at zero applied Fgd, as 
a dc voltage at the drain, while the source is grounded. 
The responsivity i?^ is extracted from the measured Au 
by using the relation: Rj^ = AuSt/iPtSa)^ where ^t is the 
radiation beam spot area and Sg, is the active area [25 . 
This definition of Ri, assumes that the whole power inci- 
dent on the antenna is effectively coupled to the graphene 
FET. In our case, the relatively high device impedance 
results in a mismatch with the antenna output (~ 100 Q 
or lower). It is thus likely that a considerable fraction 
of the radiation field is not funneled into the device. 
Therefore our R^, values have to be considered as lower 
limit estimations. Given our radiation beam spot diam- 
eter d = 4: mm, we get ^t = ttcP/A = 12.6 x 10~^ m^. 
Since the total area of our device, including antenna and 
contact pads, is smaller than the diffraction-limited area 
Sx = A^/4, we take the active area to be Sx. The in- 
duced Au was then extracted from the measured lock-in 
signal £, without any correction from the amplifier input 
load 25 : 
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FIG. 3: (Color online) (a,b) Responsivity values measured 
as a function of the gate voltage Vg at room temperature, 
for single-layer graphene (a) as function of the angle between 
the polarization axis of the beam and the antenna axis, and 
for bilayer graphene (b). The data are extracted from the 
measured photo-induced voltages across source and drain as 
described in the text. The incoming 0.3 THz beam is mod- 
ulated at 500 Hz and impinges from the free space onto the 
graphene devices mounted on the dual-in-line package with 
an incident optical power Pt = 2.1 mW. 



where the factor 2 is due to the peak-to-peak magnitude, 
the factor a/2 originates from the lock-in amplifier rms 
amplitude, and 7r/4 is the fundamental sine wave Fourier 
component of the square wave produced by the chopper. 

Fig. IsFa) plots Pjy of the SLG-based FET as a function 
of Vg, while modulating the THz source at 500 Hz. Dif- 
ferent curves correspond to different relative orientations 
between the source electric-field polarization and the an- 
tenna axis. The photoresponse drops rapidly with angle 
until it becomes almost zero when the incoming beam 
polarization is orthogonal to the antenna axis, confirm- 
ing the efficacy of our dipole antenna. The dependence of 
Au on Vg is in qualitative agreement with the behavior 
shown in Fig. [2|b) , thereby proving that in the present 
conditions our detectors are operating in the broadband 
overdamped regime. In particular, we note that the sign 
of the photovoltage changes quite abruptly in the vicinity 
of the CNP, following the switch of sign of the derivative 
da/dVg. 

We also note that an additional sign switch appears 
around Fg = in all curves. This suggests a con- 
tribution to the photoresponse of thermoelectric ori- 
gin [19l |20], arising from the presence of the ungated 
p-doped graphene regions, and subsequent formation of 
p-p-p or p-n-p junctions, depending on the gate bias 
value. It is likely that, owing to the presence of the 
antenna, the electronic distribution near the junction be- 
tween gate and source is more effectively "heated" by the 
incoming radiation, either through free-carrier (Drude- 
like) absorption or indirectly from the lattice. The re- 
sulting thermoelectric voltage should have the same func- 
tional dependence as in Eq. ^ but opposite sign, plus a 
roughly gate-independent value, determined by the ther- 
moelectric Seeback coefficient of the ungated graphene 
region |T9l |20] . It should thus change sign near zero gate 
bias, following the sign of the difference between the ther- 
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FIG. 4: (Color online) (a,b) Noise equivalent power (NEP) 
as a function of the gate voltage measured for single-layer 
graphene (a) and bilayer graphene (b) assuming that the noise 
level is dom inated b y the thermal Johnson-Nyquist contribu- 
tion N = v^4A:BT/cr. 



moelectric coefficients across the junction, as in Ref. [20] . 

Our device reaches i?j, ~ 100 mV/W, when the source 
is modulated at 1 kHz. This confirms the presence of 
competing detection mechanisms, each characterized by 
distinctive response timescales. A quantitative evalua- 
tion of each individual contribution is, however, presently 
difficult, due to their similar functional dependence. As 
inferred from Eq. Q, in the overdamped regime, R^ is 
determined by the relative change of conductivity with 
V^; in fact, in this first implementation, our R^ is almost 
one order of magnitude lower than that of nanowire FETs 
with identical antenna geometries [10 , but with an on-off 
current ratio at least one order of magnitude larger. 

Figure IsFb) plots R^ as a function of V^ in the BLG- 
based FET, while modulating the THz source at 500 Hz. 
The measured shape of the R^ curve is in excellent agree- 
ment with that predicted by the diffusive plasma-wave 
detection model [see Fig. ^d)], up to the CNP point. 
However, no change of sign is measured at the CNP, and 
a response strongly enhanced with respect to that pre- 
dicted by Eq. (fTl) appears for V^ larger than the CNP, 
when the device is in the p-n-p configuration. This sug- 
gests an additional contribution to the photovoltage, this 
time of constant sign. Its magnitude grows rapidly with 
gate bias, eventually dominating in the regime in which a 
p-n-p junction is present. We note that the photovoltage 
is noisier above the CNP, implying an additional intrinsic 
noise mechanism. We tentatively link this behaviour to 
interband transitions driven by the THz field at the p- 
n junction, with the resulting generation-recombination 
noise. Indeed, this contribution, only seen in BLG, may 
be related to the higher density of states near the Dirac 
point, or to the opening of an energy gap due to the elec- 
tric field, with the THz field possibly adding an energy 
gap modulation. Whatever the precise origin, this leads 
to Ry ~ 150 mV/W, when the source is modulated at 
1 kHz, values even larger than those resulting from the 
broadband plasma response. 

From the application perspective, the relevant figure 
of merit to characterize a photodetector is the NEP, 



which corresponds to the lowest detectable power in a 
1 Hz output bandwidth (or equivalently a 0.5 s integra- 
tion time). The noise level of FET detectors is domi- 
nated by the thermal Johnson-Nyquist contribution [9], 
expressed as A/" = ^J ^k^T j a^ and can then be extracted 
directly from the measured conductivity. The NEP is 
then simply obtained dividing this value by the detec- 
tor responsivity Ry. Figures |4ja) and (b) plot our best 
measured NEPs as functions of V^. The minimum RT 
NEP is - 200 nW/Hz^/^ for SLG, and almost one order 
of magnitude lower (- 30 nW/Hz^/^) for BLG. This dif- 
ference is expected in view of the higher responsivity of 
the latter, combined with its higher conductivity (trans- 
lating into lower noise). We note that such numbers are 
in fact upper limits, and since they refer to the signal 
power incident on the device as a whole (optical NEP), 
i.e. are not corrected for the coupling efficiency of the 
radiation into the nano-sized transistor element, which is 
likely very low in our present configuration. Care must 
be taken when comparing optical NEP values with those 
derived from the signal power actually absorbed in the 
detector (electrical NEP), used for example to describe 
the mid-IR cryogenic bolometers of Ref. [23^. Previous 
RT FET detectors based on InAs nanowires achieved 
optical NEPs - 10-^ W/Hz^/^ [10], while Si n-MOSs 
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[9 . The reason why our de- 



vices have higher NEPs stems from the higher Ry en- 
abled by the current on-off ratio of the highly mature Si 
CMOS technology, exceeding 10^. We thus expect signif- 
icant performance improvements to be possible for our 
graphene FET detectors, even in the overdamped low- 
frequency regime, for example by reducing the device 
resistance with shorter channels and raising the mobil- 
ities, improving at the same time the coupling with the 
antenna. We also remark the possibility recently demon- 
strated to reach on-off current ratios as high as 10^ in 
CVD-grown BLG FETs [26]. 

In any case, our devices are already exploitable, for 
instance, for large area fast imaging. As test object we 
used coffee capsules (aluminum foil) inside a cardboard 
box. Figure [5] shows the THz image profile, consisting 
of 200 X 550 scanned points, collected by raster-scanning 
the object in the beam focus, with an integration time of 
20 ms/point. The coffee capsules as well as the air gaps 
between them appear clearly in the THz image with a 
reasonably good spatial resolution. A transmission image 
of a fresh, green leaf is also shown [see panel (d) of Fig.ls] . 

Our SLG and BLG-based FET detectors can also be 
used to investigate more fundamental physics. For ex- 
ample, in the weakly-damped limit, they may allow to 
probe the hydrodynamic behavior of "chiral" electron 
plasmas and their non-linear instabilities [HJIS^. More- 
over, BLG FETs with narrow top gates are particularly 
suited to unravel the interplay between photo-excited 





FIG. 5: (Color online) (a-c) Photograph of a coffee-capsule 
box (a-b) and 0.3 GHz transmission mode image of it (b) 
measured at room temperature with a bilayer-graphene-based 
FET at 14 = 3 V, mounted on a XY translation stage having a 
spatial resolution of 0.5 /xm. The THz image consists of 200 x 
550 scanned points, collected by raster-scanning the object in 
the beam focus, with an integration time of 20 ms/point; (d) 
0.3 GHz transmission mode image of a leaf measured under 
the same experimental conditions. 



plasma-wave transport and intriguing phenomena such 
as chirality-assisted electronic cloaking [28] and Zener- 
tunneling-induced negative differential conductivity [29^. 
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Methods 
Fabrication and characterization procedures 

Monolayer (SLG) and bilayer (BLG) graphene flakes 
were fabricated by mechanical exfoliation from Kish 
graphite on a highly intrinsic {p = 10 kl^ cm) Si substrate 
with a 300 nm thermally grown Si02 surface layer. SLG 
and BLG flakes were identified by optical microscopy 
and by Raman spectroscopy experiments. Spatially- 
dependent Raman measurements were carried out ex- 
ploiting the 488 nm line of an Argon- ion laser at 1.5 mW. 
The laser was focused on the sample by a 0.7- numerical- 
aperture lens, giving a sub-/i7Ti illumination spot. 

Both samples were then spin-coated with e-beam sensi- 
tive resist and contact patterns were exposed by electron 



beam lithography to realize the drain and source con- 
tact pads. Gr(5 nm)/Au(80 nm) Ohmic contacts were 
then thermally evaporated onto the samples, and lift-off 
was made in heated acetone. Gontact resistivities on the 
order of 10~^ Vt cm , were measured. In a successive fab- 
rication step we defined the top gates by e-beam lithog- 
raphy, after preventively depositing via atomic layer de- 
position (ALD) a 35 nm thick Hf02 dielectric gate at 
a temperature < 200 °G. Due to the hydrophobic na- 
ture of graphene, a functionalization layer (FL) needs 
to be deposited before the ALD process. This sticks to 
graphene and provides the catalytic sites required to start 
the ALD deposition. A metallic (Gr/Au) top-gate chan- 
nel was then patterned in the middle of the graphene 
flake. We extracted a top-gate capacitance per unit area 
of about 300 — 500 nF cm , assuming a dielectric con- 
stant /^ = 13 — 19 nF cm~ of the amorphous Hf02. Fi- 
nally Raman spectra measured both before and after the 
Hf02 deposition did not show any evidence of defect- 
induced D peaks indicating no significant degradation of 
the graphene properties after oxide deposition. The elec- 
trical characterization was performed after having kept 
the devices under vacuum for 12 hours. The drain con- 
tact was connected to a current amplifier converting the 
current into a voltage signal with an amplification factor 
of 10"^ V/A. 



Diffusive transport model 

When the frequency of the impinging radiation uj sat- 
isfies the condition ujTqq <C 1, 1/tqq being the electron- 
electron scattering rate, the electron gas in the channel 
follows the laws of hydrodynamics [Tl] [27]. Electrons 
undergo many collisions among each other during one 
period of oscillation of the external field, thus establish- 
ing local thermodynamic equilibrium. The time scale 
associated with electron-electron collisions in graphene 
is on the order of ~ 1 — 10 fs [see, for example, M. 
Polini et al, Phys. Rev. B 77, 081411 (2008)]. In the 
present experiment uj/(2it) = 0.3 THz, meaning that 
cjTee ~ 2 X 10~^ — 2 X 10~^. In this regime, electrons 
move collectively establishing plasma waves in the chan- 
nel, which are however strongly damped since ujTtr <C 1. 
Here rtr is the transport time, which can be estimated 
from the field-effect mobility to be in the tens of fs. In the 
hydrodynamic regime, the photoresponse can be calcu- 
lated by using a one-dimensional model based on Ohm's 
law, 



j{x,t) = (jE{x,t) 
and the continuity equation 
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The total carrier density —en{x,t) is modulated by the 
gate voltage Vg{x,t): 



enix^t) = CVg{x^t) 



(5) 



where C is the gate-to-channel capacitance per unit area. 
The d.c. conductivity a depends on the total carrier 
density and thus on the gate voltage Vg{x^t) by virtue 
of Eq. ([5|. The system of coupled equations (|3|-([5| can 
be solved with suitable boundary conditions, which are 
imposed at the source {x = 0) and at the drain (x = Lg), 
with Lg ^ CgTtr. Here Cg is the acoustic-plasmon group 
velocity. We now seek solutions of Eqs. (|3|-([5]) of the 
form Vg{x^t) = Uq -\- Ui{x,t) + U2{x) with the following 
boundary conditions: Ui{x = 0,t) = UaCOs{ujt)^ Ui{x = 
Lg,t) = 0, U2{x = 0) = 0, and U2{x = Lg) = const. This 
leads to the following source-drain photovoltage: 



Au-- 



U^ 



da{V) 



4 a{Uo) dV 



(6) 
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